Cross-correlation Tomography: Measuring Dark Energy Evolution with Weak Lensing 
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A cross-correlation technique of lensing tomography is presented to measure the evolution of 
dark energy in the universe. The variation of the weak lensing shear with redshift around massive 
foreground objects like bright galaxies and clusters depends solely on the angular diameter distances. 
Use of the massive foreground halos allow us to compare relatively high, linear shear amplitudes 
in the same part of the sky, thus largely eliminating the dominant source of systematic error in 
cosmological weak lensing measurements. The statistic we use does not rely on knowledge of the 
foreground mass distribution and is only shot-noise limited. We estimate the constraints that deep 
lensing surveys with photometric redshifts can provide on the dark energy density f2do , the equation 
of state parameter w and its redshift derivative w'. The marginalized accuracies on w and w' are: 
a(w) ~ 0.02 f~^J 2 and a(w') ~ 0.05/ s ^^ 2 , where / s k y is the fraction of sky covered by the survey 
and a(Qdc) = 0.03 is assumed in the marginalization. Combining our cross-correlation method with 
standard lensing tomography, which has complementary degeneracies, will allow for measurement 
of the dark energy parameters with significantly better accuracy. 



I. INTRODUCTION 

Gravitational lensing provides us with the most direct 
method for probing the distribution of matter in the uni- 
verse 0. Lensing leads to a shear distortion of back- 
ground galaxy images, or a change in the surface num- 
ber density of background galaxies due to magnification. 
The measurement of the mass distribution in clusters of 
galaxies using the lensing shear and magnification 

|dl LJ. Li [Z| are now well established techniques. On larger 
scales detections of the cosmic shear signal 8j show that 
the cosmological matter distribution can also be probed 
this way. 

The variation of the lensing signal for background 
galaxies at different redshifts probes the projected lens- 
ing mass with different redshift weights in a way that 
depends on cosmology 0, E3, EH E3 Hu [II [H El 
has developed techniques for using the shear power spec- 
trum for background galaxies with photometric redshift 
information to constrain cosmological parameters, in par- 
ticular the nature and evolution of dark energy. His and 
other recent studies [H El El El El HJ have fore- 
cast the accuracy with which these parameters can be 
obtained from future weak lensing surveys, while pH |2^ | 
have developed methods to use tomography for 3-D mass 
reconstruction. Observationally weak lensing tomogra- 
phy has been applied to a galaxy cluster [23j , but further 
progress awaits multi-color imaging surveys that can ob- 
tain photometric redshifts of background galaxies. 

The lensing shear (or magnification bias) can be 
used to cross-correlate large foreground galaxies (asso- 
ciated with the lensing mass) with background galax- 
ies which are lensed. In this paper we will use such 
cross-correlations as an alternative way of doing lens- 
ing tomography. We use a particularly simple cross- 
correlation statistic: the average tangential shear around 
massive foreground halos associated with galaxy groups 
and galaxy clusters. We show that cross-correlation to- 



mography measures ratios of angular diameter distances 
over a range of redshifts. The distances are given by inte- 
grals of the expansion rate, which in turn depends on the 
equation of state of the dark energy. Thus it can be used 
to constrain the evolution of dark energy. We estimate 
the accuracy with which dark energy parameters can be 
measured from future lensing surveys. 



II. FORMALISM 



We work with the metric 



ds 2 



a 2 [-(1 + 2<j))dT 2 + (1 - 2<j)){d X 2 + rW)] , 



(1) 

where we have used the comoving coordinate x an d 
a(r) = (1 + z)^ 1 is the scale factor as a function of 
of conformal time r. We adopt units such that c = 1. 
The comoving angular diameter distance r(x) depends 
on the curvature: we assume a spatially flat universe so 
that r(x) = X- The density parameter f2 has contribu- 
tions from mass density Q m or dark energy density Qdc, 
so that Q — Q m + Qdc- The dark energy has equation of 
state p = wp, with w = — 1 corresponding to a cosmolog- 
ical constant. The Hubble parameter H(a) is given by 



H{a) = H 



n m a~ 3 + n dc e- 3 t? dlna ' {1+w{a,)) 



(2) 



where H is the Hubble parameter today. The comoving 
distance x(a) is 



X(a) 



da' 



! a' 2 H(a') : 



(3) 



The lensing convergence is given by the weighted projec- 
tion of the mass density 

k{4>) = -n m / d x g{x) , (4) 
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where the radial weight function gb(x) can be expressed 
in terms of r(x) an d the normalized distribution of back- 
ground galaxies Wb{x) 



9b(x) = r (x) 



r{x' - X) 
r(x') 



W h (x'W , 



(5) 



where xo is the distance to the horizon. For a delta- 
function distribution of background galaxies at Xbj this 
reduces to g h (x) = KxMXb - xVKXb)- 

We consider the lensing induced cross-correlation be- 
tween massive foreground halos, which are traced by 
galaxies, and the tangential shear with respect to the halo 

center (denoted 7 in this paper): uj x (9) = (^Snf($)j($')J 

where ri{(<p) is the number density of foreground galax- 
ies with mean redshift (z { ), observed in the direction <j> in 
the sky and Sri{((f>) = (n{(<j)) — n{)/n{. The angle between 
directions 4> and 4>' is 9. The cross-correlation is given by 
MM: 



uj x (9) = 67r 2 rj K 



a(x) 



< / dkkP hm {x,k) J fi [kr(x)9] 
h 



(6) 



where Phm(Xi^) is the halo-mass cross-power spectrum, 
and Wf is the foreground halo redshift distribution. The 
Bessel function J M has subscript // = 2 for the tangential 
shear and /1 = for the convergence (from the relation 
7 (0) = -1/2 dR(9)/d\n9, 25]). The measurement of 
the mean tangential shear around foreground galaxies is 
called galaxy-galaxy lensing. We will consider a general- 
ization of this to massive halos that span galaxy groups 
and clusters. 

If the foreground sample has a narrow redshift distri- 
bution centered at x = Xit then we can take Wf to be a 
Dirac-delta function and evaluate the integral over x- All 
terms except gt(xf) are then functions of Xf , the redshift 
of the lensing mass. The coupling of the foreground and 
background distributions is contained solely in gb(Xf)- 
Hence if we take the ratio of the cross-correlation for two 
background populations with mean redshifts Z\ and Z2, 
we get 



wi(9) = gi(xf) 
w 2 (0) 92 (Xf)' 



(7) 



where w\ , g\ denote the values of the functions for the 
background population with mean redshift Z\. In the 
limit that the background galaxies also have a delta- 
function distribution, this is simply a ratio of distances 



wi{9) = r(xi - Xf)A(xi) 
w 2 {6) r(x2 - Xf)/r(x2) 



(8) 



The above equations show that the change in the cross- 
correlation with background redshift does not depend on 



the galaxy-mass power spectrum, nor on 9. We can sim- 
ply use measurements over a range of 9 to estimate the 
distance ratio of equation (jHJl for each pair of foreground- 
background redshifts. The distance ratio in turn depends 
on the cosmological parameters fide, w, and its evolution 
w'. We evaluate next how the cross-correlation defined 
above can be used to constrain the dark energy parame- 
ters. This approach is complementary to standard lens- 
ing tomography, where the variation of power spectrum 
of the shear with background redshifts is used to mea- 
sure an integral of distance and growth factors over red- 
shift. With the cross-correlations, there is no dependence 
on growth factor, but a precise dependence on distances 
(since they are not integrated over redshift). 



III. RESULTS 
A. Signal-to-noise estimate 

A simple way to estimate the signal-to-noise for the 
cross-correlation approach is to regard the foreground 
galaxies as providing a template for the shear fields of 
the background galaxies (G. Bernstein, private commu- 
nication). For a perfect template (i.e. for high density 
of foreground galaxies and no biasing), the errors are 
solely due to the finite intrinsic ellipticities of background 
galaxies. Thus the fractional error in our measurement 
of the background shears is simply 



1 



1' 



(9) 



The total number of background galaxies is Atotai = 
n g A = n g /sky^4sky, where A is the survey area, the lens- 
ing induced rms shear (7) rm s — 0.04 0, and the intrinsic 
ellipticity dispersion a e — 0.3. This gives 



1 



0.2 x 10" 



100 

n a 



1/2 



0.1 

/sky 



1/2 



(10) 



where the number density n g has units per square ar- 
cminute. Thus for the fiducial parameters / s k y = 0.1 and 
n g = 100, one can expect to measure the background 
shear to 0.1% accuracy at about 5-cr. We find that such 
a signal corresponds to changes in w of a few percent; 
hence this is the approximate sensitivity we expect in 
the absence of systematic errors. 



B. Models 

To construct a simple and observationally robust cross- 
correlation statistic, we restrict the foreground sample to 
galaxy clusters and large galaxy groups. With deep imag- 
ing surveys, most of these massive objects are expected to 
be identified out to z ~ 1. We use the tangential compo- 
nent of the shear for two background samples at different 
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redshifts, inside apertures of size # ap ~ 3', as our estima- 
tor of the distance ratio of equation (JJJ . The signal from 
all massive halos in a given redshift bin will be stacked for 
each background sample. The tangential shear around a 
foreground halo is (7(6*)) = (E(0) - £(0))/£ crit where 
S C rit is the geometrical factor that depends on angu- 
lar diameter distances. To compute the average shear 
around halos within a mass range, one replaces the pro- 
jected mass density S by the projection of the halo-mass 
correlation function, so that the above equation is equiv- 
alent to equation © in the limit of thin redshift distribu- 
tions. Wc will assume that the massive foreground halos 
at z < 1 have spectroscopic or accurate photometric red- 
shifts, since they generally contain at least one bright 
galaxy. Each background galaxy sample is described by 
a redshift distribution determined using photometric red- 
shifts (which can have relatively large statistical errors). 

The average tangential shear around foreground galax- 
ies and clusters has proven easier to measure than shear- 
shear correlations since it is linear in the shear and has 
a distinct signal that can be distinguished from system- 
atic contributions (for a wide field survey it is unlikely 
for a systematic to produce an averaged tangential shear 
around foreground galaxies). Moreover, in comparing the 
signal for two background populations, it is a great ad- 
vantage that the shear amplitudes are compared for the 
same set of apertures on the sky: the relative amplitudes 
are insensitive to spatial variations of the point spread 
function, which are the dominant systematic for shear- 
shear correlations. 

To compute the accuracy with which parameters can 
be measured, we need to know the signal and the noise. 
The signal is given by equation JJjJ in which the halo- 
mass cross-correlation can be accurately computed using 
the halo model of large-scale structure (see |2g for a re- 
view). The model specifies n(m,z), the comoving num- 
ber density of halos with mass m, the bias parameter of 
halos, and the density profile of a halo for given mass and 
redshift. Details of the model ingredients relevant to this 
study are given in [53, ■ 

The two key quantities for each redshift bin are the 
number density of halos and the mean shear for each ha- 
los mass. We choose the mass range 4 x 10 13 < to/ Mq < 
10 15 and aperture size 9 ap = 2'— 3' corresponding roughly 
to the virial radii of halos over the dominant redshift 
range. We also set a constant inner aperture 8 — 0.1' to 
exclude the strong lensing regime and the luminous parts 
of the lensing halos. With foreground redshift slices of 
width Az = 0.1, the halos in each slice cover about 1/10 
the survey area. Thus using all the foreground slices over 
< Zf < 1, close to the entire survey area is covered. 

We take the background redshift distribution to be 

given by dN/dz ~ z a exp — (z/zo) b , with a = 2 and 

b = 1.5 and zq — 1 giving a mean redshift of 1.5. Wc 
normalize it to give a total number of density n g = 100 
per square arcminute. It is split into two background 
samples, between Zf and z = 3.0, with about half the 
galaxies in each sample. The halo model is used to ob- 



tain the mean shear for each lens-source redshift distri- 
bution by summing the contributions of halos over the 
chosen mass range. While we have used the halo model 
to compute the expected signal-to-noise, the inferred pa- 
rameter values do not rely on knowing the masses or any 
properties of the halos. 



C. Dark energy parameters 

We perform a \ 2 minimization over the mean shear 
amplitudes at the two background distributions for each 
foreground slice to fit for the dark energy parameters. 
The time dependence of w is parameterized as w = wq + 
w a (l — a) |31| . For comparison with other work we will 
compare w a to w' defined by w — u>q + w'z. Since the w' 
parameterization is unsuitable for the large redshift range 
we use, any comparison with it can be made only for a 
choice of redshift. For example, at z = 1, which is well 
probed by our method and is of interest in discriminating 
dark energy models .31], w a — 2w' . For foreground slices 
labeled by index I and two background samples by 1 and 
2, the x 2 is given by 



Y 



E 



1 - 



R°(z h zi,z 2 ) 
R(zi,zi,z 2 ) 



(11) 



where R is the distance ratio of equation Q for given 
values of fide j w and w a , and R is the fiducial model 
with fide = 0.7, w = —l,w a = 0. The weights Ui are 
given by 



u; 



2mfiA ( 7 >?i 2n 2 f t A ( 7 >f 



(12) 



where fi is the fraction of the survey area A covered 
by halo apertures in the l-th lens slice. The factor of 2 
in the denominator arises because we are using only one 
component of the measured ellipticity whereas of denotes 
the sum of the variances of both components. 

The results are shown in Figures ^ an d [21 Figure 2] 
shows the constraints in the fide — w plane. The ellipses 
show the 68% confidence region given by A\ 2 = 2.3. 
The elongated inner contour is for fixed w a — 0, while 
the outer contour marginalizes over u> a ; an external con- 
straint on fide from the CMB and other probes is as- 
sumed, corresponding to er(fid c ) = 0.03. 

Figure |3 shows the constraints in the w — w a plane 
if fide is fixed, or marginalized over with cr(fidc) = 0.03. 
The corresponding accuracy on w and w' can be scaled as 
/sky /2 - For °"( fi de) = 0.03, we obtain a(w) ~ 0.02/ s 7 1/2 



— 1/2 

and a(w a ) ~ 0.1/ sky (at z = 1 this is equivalent to 

— 1/2 

<j(w') ~ 0.05/ sky ). Note that the value of a on a param- 
eter is given by projecting the A\ 2 — 1 contour on the pa- 
rameter axis, which is smaller than the projection of the 
inner contour in Figure [3 with A\ 2 = 2.3. The scaling 
with / s ky in the parmameter errors comes from the num- 
ber of background galaxies (not sample variance), hence 



'sky 
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FIG. 1: Contours in the Qdc — w plane for the fiducial lens- 
ing survey with / s k y = 0.1. The inner contour assumes no 
evolution of dark energy, w a = 0, while the outer contour 
marginalizes over w a ; both employ a prior on fide correspond- 
ing to a(Qdc) = 0.03 (see text). The 68% confidence interval 
is shown in each of the contours. 



FIG. 2: Contours in the w — w' plane for the fiducial lensing 
survey with / a k y = 0.1, as in Figure Q The inner contour 
assumes Q^c ~ 0.7, while the outer contour marginalizes over 
fide with the prior cr(fidc) = 0.03. Note that the parameter 
w a = 2 w' at z = 1 (see discussion in the text). 



IV. DISCUSSION 



for given / s k y varying the depth of the survey scales the 

— 1/2 

errors roughly as n g . The results we have shown are 
for the fiducial redshift z = 0. A different choice of the 
fiducial redshift changes the relative accuracy on w and 
w a somewhat, because the degeneracy direction in the 
three parameters changes. A detailed exploration of dif- 
ferent models of w(a) with finer bins in the background 
redshift distribution would be of interest. 

The parameter constraints obtained above are consis- 
tent with the crude signal-to-noise estimate made above 
in equation (|10|l . This can be seen by using representative 
numbers for the halos used in the cross correlation. For 
m = 10 14 Mq, the number density n(m) ~ 0.01/arcmin 2 
per In interval in m integrated out to Z = 1, while the 
mean shear inside 6* ap ~ 2' is (7) ~ 5% for Zf — 0.3 and 
Zh = 1 28] . Thus using halos with a mass range of over 
one decade centered on m = 10 14 M Q gives close to full 
sky coverage. For relevant scales, the shear profile around 
the halo center is close to the isothermal form: 7 oc 1/8. 
The signal-to-noise is then constant per log 8 since the 
noise variance scales as the area (oc 8 2 ). We will use the 
simple scheme of stacking all halos, and averaging over 
the mean shear inside apertures, to estimate the signal. 
The noise per unit area is a t / \J2n g . Using the numbers 
given above for (7) and the noise gives a signal-to-noise 
estimate in agreement with equation I|1U|I as well as the 
results shown in Figures ^ an d 121 



We have presented a method of lensing tomography 
that uses foreground halos associated with the lensing 
mass to measure angular diameter distances over a range 
of redshifts, < z < 3. It offers an observationally robust 
probe of dark energy and its evolution over this redshift 
range. The accuracy achievable with a lensing survey 
that covers 1/10 of the sky (with parameters close to 
that of the survey proposed with the LSST telescope) is 
potentially one of the most promising for proposed sur- 
veys in the coming decade. It is complementary with 
constraints from Type la supernovae in that the redshift 
coverage is broader while the distance factors probed 
are similar [2(|. Combining cross-correlation tomogra- 
phy with the standard shear power spectrum tomography 
improves the constraints significantly, since the contour 
ellipses are oriented differently 0] - this is especially 
appealing because the two techinques use the data in dif- 
ferent regimes, the linear and strongly nonlinear regime. 
Further it would allow the constraints on the dark mat- 
ter, such as on neutrino mass, from shear power spectrum 
tomography to be improved by providing independent in- 
formation on the geometry. 

An attractive feature of cross-correlation tomography 
is its observational feasibility. The main practical ad- 
vantages are: (i) large shear values (~ 1 — 10%) around 
massive halos are used, and (ii) linear shear amplitudes 
are compared in the same apertures on the sky. The 
statistic is thus largely insensitive to the primary limit- 
ing systematic error in weak lensing surveys: the varia- 
tion of the point spread function over the field view (to 
first order the psf variation averages to zero for a linear 
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statistic) . The imaging requirements are far less stringent 
than for other cosmological applications of lensing, such 
as standard tomography, which rely on the measurement 
of quadratic shear correlations at large angles. The on- 
going CFH Legacy survey should allow for testing of the 
method and parameter constraints that can be checked 
with limits from the CMB. The limiting systematic er- 
ror for most surveys is likely come from the photometric 
redshifts of the background galaxies. While large statis- 
tical errors do not affect the parameter constraints signif- 
icantly, systematic errors could bias the inferred distance 
ratios. Calibrating a fair sub-sample of the photomet- 
ric redshifts with spectroscopic redshifts, and testing for 
effects such as possible correlations between measured el- 
lipticites and the sizes or surface brightness of galaxies, 
should allow us to safeguard against such biases p9|. 

We have assumed in this study that the intrinsic ellip- 
ticities of background galaxies dominate the statistical 
errors. Other sources of errors are projection effects due 
to matter at different redshifts from the foreground halo 
and intrinsic correlations in the ellipticies of background 
galaxies. Our choice of averaging the tangential shears 
around halos should cause these errors to not bias the 
result, but a careful study is needed to quantify the con- 
tribution to the statistical error. A possible bias can be 
introduced if the foreground redshifts have large errors 
and the halo selection, based on visible galaxies, is per- 
versely correlated with dark energy evolution. Statistical 
accuracy in the foreground redshifts of better than 0.1 in 
z should rule out such a bias due to cosmological evolu- 



tion, but again a more detailed study is needed. 

Our implementation of halo-shear cross-correlations is 
not optimal in that we have used only a fraction of the 
measured shapes for each lens slice, and only two bins of 
the background galaxies (motivated in part by the finding 
that this provides most of the information in shear power 
spectrum tomography [T^| ) . The parameter accuracy can 
therefore be improved with an optimal scheme in which 
we use finer binning or the actual redshifts (or photomet- 
ric redshift probability distributions) of foreground and 
background galaxies rather than binning them. It will 
also be of interest to study different models of w(z) rather 
than the monotonic parameterization used here, or to ob- 
tain model independent constraints on the expansion his- 
tory directly [30l l3l| . A joint analysis with shear-shear 
correlations is needed to quantify the improved precision 
one can expect. We have restricted ourselves to the weak 
lensing shear; additional information can be obtained by 
using the strong lensing signal expected in some fraction 
of the galaxy and cluster halos (for related studies with 
lensing arcs from galaxy clusters, see [32, 03) an d 
using the magnification signal. 
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